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Easy steps through essential Hamiltonian dynamics are outlined, from 
necessary definitions of Poisson brackets and canonical transformations, to 
a quick proof that Hamiltonian evolution is made by canonical transforma- 
tions, the quickest proof of Liouville's theorem, and on to Poincare-Cartan 
integral invariants and completely integrable dynamics, making room, pro- 
viding tools, and setting the stage for more recent developments. 
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I. Introduction 



At a few key points, my teaching of classical mechanics relies on calcu- 
lations that I do differently from standard texts. Strung together with bits 
of familiar material, they form an alternative route through Hamiltonian dy- 
namics that I find advantageous. It fits more easily into a course that moves 
from traditional to chaotic dynamics. This requires that I provide notes for 
students. Here I describe points that may be helpful for teachers as well. 
They supplement what I have found in the standard texts. I take care to 
refer to the texts I have found most helpful and those that are closest to 
what I do. 

New developments in classical mechanics have given Hamiltonian 
dynamics a new role in our teaching. It helps prepare us to work with the 
methods and understand and appreciate the results of chaotic dynamics. 1-4 
To make room for new developments, we have to delete or reduce some 
traditional topics. Hamiltonian dynamics involves an accumulated wealth of 
mathematical development that can be shortened without sacrificing either 
the basic concepts or the results that illuminate contemporary work. We 
do not need to know much about canonical transformations, for example, 
or anything about their generating functions, to work with the canonical 
transformations of Hamiltonian evolution. Using the way they depend on 
time can be enough. 

We can edit Hamiltonian dynamics for economy of expression and 
aim it in directions of current interest. We can, if we choose, step lightly 
over some traditional theory and save time for solving problems, or have 
time to introduce modern geometrical views. 1 Minimal steps through essen- 
tial Hamiltonian dynamics are outlined here, from necessary definitions of 
Poisson brackets and canonical transformations, to a quick proof that Hamil- 
tonian evolution is made by canonical transformations, the quickest proof of 
Liouville's theorem, and on toward contemporary goals, to Poincare-Cartan 
integral invariants, providing tools that can be used with surface-of-section 
maps, 5 and to completely integrable dynamics, setting the stage for use of 
action-angle variables and discussion of the effect of perturbations on reso- 
nant and non-resonant tori and the KAM theorem. 6 Generating functions 
are particularly useful for proving that the transformation to action-angle 
variables is canonical, but that is proved here simply by calculating Poisson 
brackets. We do not need a power tool for this one job; we can do it by hand. 
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II. Poisson Brackets 



We work with canonical coordinates qi, q^ ...q^ and momenta pi,P2, ■ ■■Pn- 
The Poisson bracket is defined by 



n=l ^ 



/ dF dG dF dG 



dq n dp n dp n dq n 



(2.1) 



for any two functions F and G of the canonical coordinates and momenta. 
It is linear for both F and G . It is antisymmetric; [G, F] is — [F, G] . These 
properties are obvious. It is an easy exercise to prove that the Posisson 
bracket also satisfies the Jacobi identity 

[[F, G},H] + [[G, H],F] + [[H, F},G} = 0. (2.2) 

The canonical coordinates and momenta themselves have Poisson brackets 

[q m ,q n } = o , [p m ,p n ] = o 

[q m ,Pn]=5mn (2.3) 

for m,n = 1,2, ... , N. Hamilton's equations 

dq n dH dp n dH 
dt dp n ' dt dq n 

can be written as 



(2.4) 



Indeed, Hamilton's equations and the definition of the Poisson bracket imply 
that for any function F of the canonical coordinates and momenta 

*±(fMf"L) = lF , m . (2 . 6) 

The Hamiltonian H is a function of the canonical coordinates and momenta 
that may be different for different t; it may be a function of the q n ,p n and t. 
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III. Canonical Transformations 



The 2N-dimensional space of points specified by the canonical coordi- 
nates and momenta is called phase space. We consider changes of coordinates 
in phase space. Let Q n ,P n , for n = 1,2,. ..,N, be functions of the qk,Pk that 
determine the qk,Pk as functions of the Q n ,P n , so the Q n ,P n label the points 
of phase space as well as the qk,Pk ■ We say that the Q n ,Pn are canonical co- 
ordinates and momenta, and the transformation from the qk,Pk to the Q n ,P n 
is a canonical transformation, if 

[Qm; Qn\ , \Pmi Pn] 

[Qmi Pn\ = $mn (3-1) 

for m,n = 1,2,. ..N. Then the Q n ,Pn can be used in place of the q n ,p n to 
calculate Poisson brackets; for any functions F and G of the phase-space 
variables, 

_ dF n dG dF pi dG 
1 ' 1 ~ ^ { dQ^ [Qm > Qnl Wn + dQ^ [Qm > Pnl dP n 



TL-rn 



I 1 mi ^cn\ an L m > n 



dP m ^—'dQ n 1 dP m l ~ m) nl dP, 



n 



f)D f)P r)P PHI I ' 1 ' ' 



dQ n dP n dP n dQ n 
From (2.6) and (3.2) we see that 



dt 1 dP, 



n 



Hamilton's equations are the same for the Q ni P n as for the q n ,p n - The Hamil- 
tonian is the same thing written differently as a function of different variables. 
Conversely, if this is true for any Hamiltonian, the transformation from the 
qt,Pk to the Q n ,Pn must be canonical; with (2.6) used to write the time 
derivatives as Poisson brackets, Hamilton's equations for the Q n ,P n give 

fill PjZJ 

v).,m = m . [ft.fli = -«j-. (3.4) 

which imply (3.1) in the particular cases where H is Q m or P m . 
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IV. Hamiltonian Evolution 



When the dynamics is described by Hamilton's equations, the evolu- 
tion in time is made by canonical transformations. Let Q n and P n be the 
functions of the qk,Pk and t that are the solutions of the equations of motion 

^f = [Qn,H] , ^ = [Pn,H] (4.1) 

specified by the boundary conditions 

Q n (t = 0) = q n , P n (t = 0)=p n . (4.2) 

The canonical coordinates and momenta q n and p n at time zero evolve to 
Q n ,P n at time t. If their sets of values are q and p at time zero, their values 
at time t are Q n (q, p, t) and P n (q, p, t). 

Explicitly, when the Hamiltonian does not depend on time and the series 
converges, 

Qn = Qn + t [q n , H] + -t 2 [[q n , H] , H] 

- ■ ... 7,..// ....// • ... (4.3) 

the bracket with H being taken k times in the term with t k ; this and the 
corresponding formula for P n do satisfy the equations of motion (4.1) and 
the boundary conditions (4.2). 

The Q n ,Pn are canonical coordinates and momenta. The transformation 
from the q n ,p n to the Q n ,P n is a canonical transformation. We prove this by 
showing that the Q n ,P n satisfy the bracket relations (3.1). Using the Jacobi 
identity (2.2) , we get 

J t [Qm, Pn] = [[Qm, H] , P n ] + [Q m , [P n , H}} 

= [[Q m ,P n ],H}. (4.4) 

The Poisson bracket [Q m , P n ] must be the function of the qk,Pk and t that is 
the solution of this equation of motion specified by the boundary condition 

[Q m , Pn] (t = 0) = [q m , q n ] = 5 mn , (4.5) 

which is [Qm, P n ] = 8mn- Explicitly, for a time-independent Hamiltonian, 

[Qm, Pn] = [q m ,Pn] + t [[q m , Pn] , H] + ••• = S mn . (4.6) 
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The other bracket relations (3.1) can be proved similarly. 

The notation used here distinguishes the coordinates and momenta 
Q n ,Pn at time t from the q n ,Pn at time zero. We will use this notation 
again when it is helpful, to consider the coordinates and momenta at time 
t as functions of the initial values. More often, since now we know that the 
coordinates and momenta are equally canonical at each time, we will follow 
the usual practice and work with canonical coordinates and momenta that 
depend on time, as we did in writing Hamilton's equations in Section II; this 
will be indicated by absence of Q n , P n . 



V. Integral Invariants 

A point in phase space marks a set of values for the canonical coordinates 
and momenta. It moves through phase space as they change in time. The 
character of this motion reflects the structure of the dynamics. Consider a set 
of points that occupy a volume in phase space. If the dynamics is described 
by Hamilton's equations, the size of the volume the points occupy does not 
change in time as the points move, even though its shape generally does. 
This property of Hamiltonian dynamics is called Liouville's theorem. 

The proof is very simple. 8 At any time, the time derivative of the volume is 
the integral over the surface surrounding the volume of the normal component 
of the velocity of the motion of the points on the surface. That is the integral 
over the volume of the divergence of the velocity, which is zero, because for 
Hamiltonian dynamics the divergence is 



A / _d_dqn _d_dPn \ = y-v ( d 2 H d 2 H i 

2-^1 \ Fin rli Ft-n rli I 2^, \ Fin Fir, Fir, Fin ' ^ ' ' 



. x dq n dt dp n dt J ^ \dq n dp n dp n dq n 

We see that this holds for a time-dependent Hamiltonian; the Hamiltonian 
can be a different function of the qk,Pk at each time. 

Finding quantities that do not change in time helps build a picture of 
the dynamics. The volume of Liouville's theorem is one. Others can be 
constructed as integrals along curves. 

Consider a curve of points in phase space. Let be a parameter that 
measures location along the curve. The phase-space coordinates q n ,p n of 
points on the curve are functions of 6. As a point on the curve moves in time, 
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its phase-space coordinates change from q n (9),p n (9) at time zero to 



at time t, where Q n and P n are the solutions of the equations of motion 
(4.1) that satisfy the boundary conditions (4.2). We write q for the set of qu 
and p for the set of pu- The phase-space coordinates of the moving point 
are functions of 9 and t. We can let points at different 9 move for different 
intervals of time t. If we make the time interval zero at the ends of a curve 
segment, the moved segment and the original segment will form a closed loop, 
as shown in Figure 1. We can replace the original segment with another 
moved curve to make a closed loop on which the phase-space coordinates 
change as functions of both 9 and t all the way around. In this way we can 
generate any closed loop in phase space that surrounds an area of surface of 
tracks of motion in time, an area of surface swept by a curve segment moving 
in time. 

For any such loop, 



Q„(q(0),p(0),i) 



P„(q(0),p(0),t) 



(5.2) 




n=l 



(5.3) 



the integral being around the closed loop. We can see this from 





(5.4) 



with 





(5.5) 



n=l 



n=l 



the integrals involving A being around the closed loop in the 9, t plane and 
over the area it encloses; because for Hamiltonian dynamics 




JV 



N 



n=l 



n=l 
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N 



dH dQ 



N 



dP n dH 
89 dP n 




dQ n 00 



n 



)-E 



(5.6) 



n=l 



n=l 



This implies that when a closed loop in phase space is changed by letting 
the points on it move for different intervals of time, the integral around the 
loop 



is not changed. We can see this by letting our original curve be a closed 
loop and letting the points on it move for different intervals of time to form 
a moved loop. Connecting the two loops along two closely adjacent tracks 
of motion, as in Figure 2, yields a closed loop to which (5.3) applies. It 
surrounds a surface of tracks of motion, a surface swept by a curve moving in 
time. Since the integrals in opposite directions along the two tracks of motion 
cancel, the integrals in opposite directions around the two loops must cancel, 
which means the integrals in the same direction around the two loops are the 
same. 

These results are important for surface-of-section maps. 5 The quanti- 
ties (5.7) that do not change in time are called Poincare-Cartan integral 
invariants. 9 When the Hdt term is absent, they are called Poincare integral 
invariants. 10 The Hdt term will be absent if there are no time differentials or 
if H can be taken outside the integral. We will use (5.3) in the next section. 




n=l 



(5.7) 



8 



VI. Complete Integrability 



A canonical transformation can make the equations of motion so simple 
that their solutions become trivial. This is possible if the dynamics are 
completely integrable. That means there are N constants of the motion, 
which we call I n with n running from 1 to N, that in addition to 

[I n ,H] = (6.1) 

which says they are constants, satisfy 

[J ro ,/„]=0 (6.2) 

for m,n — 1,2,...N. It is assumed also that these I n are functions of the 
qki Pk that determine the pk as functions of the q n , I n in a sufficiently relevant 
region of phase space, so the q n , I n label those points of phase space as well 
as the q k ,Pk- 

These assumptions of complete integrability imply, as we will show, that 
there are Q n for which the Q n , I n are canonical coordinates and momenta. 11 
The P n are the I n . They are all constants of the motion. That means Hamil- 
ton's equations give 



dH dP r , 



dQ n dt 

Then if is a function only of the P n and 

dQ n dH 

~dT ~ dK 



0. (6.3) 



(6.4) 



is a constant of the motion. The equations of motion are solved. The P n 
are constant in time. The Q n change at constant rates; they are linear func- 
tions of time. Action-angle variables are a particular case of these canonical 
coordinates and momenta Q n and I n for completely integrable dynamics. 

We will take two steps to prove that there are canonical coordinate part- 
ners Q n for canonical momenta I n . The first step, needed to define the Q n , 
is to show that there is a function S of the q n , I n such that 

Tr-=Pn- (6.5) 
dq n 

It is an action integral. Let 

5(q,I)= / 5> fe (q',I«. (6.6) 
J i° k=i 
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We write q for the set of q n and I for the set of I n . For each fixed set of values 
of the I n , the integral gives S as a function of the q n . The integral is along 
a path on the surface in phase space where the /„ are constant, starting at 
the point where the q have values qo- To establish that S is well defined, we 
need to show that the integral is not changed by variations in the path. 

We have described the motion in phase space generated by the Hamil- 
tonian H for evolution in time. Now we consider also the motion that is 
generated the same way when the Hamiltonian is replaced by one of the I n 
or a linear combination of the I n . From the bracket relations (6.2) we see 
that this motion does not change any of the I m . It stays on the surface in 
phase space where the I m stay constant. The velocity vector for the motion 
generated by I n , the vector that points along the direction of the motion in 
phase space and measures the rate of the motion with respect to a parameter 
s along the curve of the motion, has components 

dqi _ din dq N _ dl n 

ds dpi ' ds dpjq 

dpi = _dln dm = _dln_ , _v 

ds dqi ' ds dqN 

We are assuming that the q n , I n are coordinates for a region of phase space. 
For each fixed set of values of the q n , the I n must be coordinates for an N- 
dimensional space. They must vary in N different directions. Their gradients 
must be N linearly independent vectors. From (6.7), which relates the com- 
ponents of the gradient of I n to the components of the velocity vector for the 
motion generated by I n , we see that the N velocity vectors for the motions 
generated by the I n must be linearly independent. The velocity vectors for 
the motions generated by the I n and linear combinations of the I n point in 
all directions on the surface where the I m are constant. 

Now we can prove that the integral (6.6) is not changed by variations in 
the path on the surface of constant I m . If an infinitesimal segment of the path 
is moved on the surface of constant I m , the moved segment and the original 
segment form a closed loop, as in Figure 1. Now the motion is generated 
by a linear combination of the I n instead of H. We can use (5.3) with this 
replacement. The Hdt term becomes zero, because the linear combination 
of the I n is a constant that can be taken out of the integral, and the change 
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in the parameter t or s around the closed loop is zero. Therefore 



r N 

i^P fc dQ fc = (6.8) 
fe=i 

for the closed loop, which means the integral (6.6) is the same for the new 
path as for the old. We conclude that S is well defined. It clearly satisfies 
(6.5). 

That allows us to define the Q n . Let 

dl 

Qn = (6.9) 

The final step is to show that the Q n and I n are canonical coordinates and 
momenta, that the bracket relations (3.1) are satisfied with I n for P n and 
these Q n . 

Using (6.2), (6.5) and (6.9) and taking partial derivatives of S, the Q n 
and the p n with respect to the qk,h and partial derivatives of the I n with 
respect to the qk,Pk, we get 

[Q m , 4] = d -%r fe' J «] + 7?r J ™] 



dL 

3=1 3 3=1 3 



^ d 2 S dl n ^> dpj dl n 
dqjdlm dpj dl m dpj 

= dlndpj_ dl n dqj 
d Pj dl m d qj dl m 

dl 

= W^=6mn, (6.10) 
dl m 

in this context dqj/dl m is zero because the derivatives of I n are with respect 
to the qj and pj and the derivatives of the pj and qj are with respect to the 
q n and I m . Similarly, we get 
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N dQ dQ 

[Qm, Qn] = Yl -fT 1 ft! 



^ <)<,, dq 3 

E dQrn , dlj dQ n dQ n 
FlT- ^ Fl-n' Fin- 



l dl; dpj dqj dl m 

d 2 S dli d 2 S d 2 S 
< * Pit Pit Pin - Pin Fit 



N 



< x dlidl m dpj dqjdln dl m dl n 



v Q2s dlidpj 92 s 

2^/ pit pit Th^-Trr ^ 



^ dlidlm dpj dl n dl m dl n 

+ ^ = (6.1!) 



dl n dl m dl m dl n 



by using again the last steps of (6.10) which show that in this context 

W,ar m =k "- (6 - 12) 
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Figure 1. A closed loop made by points on a curve 
segment moving for different intervals of tune . 
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Figure 2. A closed loop made by an initial loop and a 
moved loop connected along closely adjacent tracks of 
motion. 
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